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Abstract 

We calculate the masses and the strong decays of orbitally excited states Bsq, B'^i, B^i and 
Bs2 by the improved Bethe-Salpeter method. The predicted masses of Bsq and B'^^ are Mb^q = 
5.723 ±0.280 GeV, Mr/ = 5.774 it 0.330 GeV. We calculate the isospin symmetry violating decay 

si 

processes Bsq — s- Bgir and B'^^ B*7r through ir^ — rj mixing and get small widths. Considering 
^ the uncertainties of the masses, for BgQ and B'^^, we also calculate the OZI allowed decay channels: 

Bso BK and B'^^ -> B*K. For Bsi and Bs2, the OZI allowed decay channels Bgi — > B*K, 
D . Bs2 ~^ and Bs2 B*K are studied. In all the decay channels, the reduction formula, PCAC 



relation and low energy theorem are used to estimate the decay widths. We also obtain the strong 
coupling constants Gb^qBstt, G^^gBK^ Gb'^_^b*tt, ^b'^-^B'tt, G'g^jS'i?) ^b'^^B'K^ Gb^^b-k^ ^b,iB*k^ 
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I. INTRODUCTION 



The heavy-light mesons play an important role in hadron physics. During the past several 
years, many heavy-light mesons were observed in experiments. In the Particle Data Group 
(PDG) table [Ssl, there are four P-wave charm states ^^(2400)°, Di(2420)°, L>i(2430)°, 
D^(2460)°, and two P-wave bottom states Pi (5721)°, P|(5747)°. For the P-wave bottom- 
strange states, Bsi and P*2 are observed by the CDF Collaboration in 2008 [11]. Later the 
DO Collaboration also reported P*2 12|. Meanwhile DO Collaboration indicated that Bgi 
was not observed with the available data set. 

In the heavy quark effective theory (HQET) 13|], for the heavy-light meson system, the 
angular momentum of light quark ji is a good quantum number when the heavy quark have 
TTiQ — )■ oo limit. They are jf = | H doublet (0~,1~) with orbital angular momentum 

1^ T doublet (1^,2+) with orbital angular 



L = 0; jf = 1^ 5 doublet (0+, 1+) and jf = 
momentum L = 1. The DO and CDF indicated that P^i (5830) and P*2(5840) correspond 



to the states with = 1+ and = 2+ in T doublet [uJ, ll2|. While for the P*j state 
with S doublet = (O"*", l"*") do not have the experimental evidence. For the charm states, 
P'*(2400)° and , Pi (2430)° are the |+ S doublet (0+, 1+), Pi (2420)° and P*(2460)0 belong 
to 1^ T doublet (l+,2+). Pi(5721)0 P*(5747)° also belong to |^ T doublet (l+,2+). 
The observations of these P-wave mesons inspire our interest in their nature. There are 



many theoretical approaches are used to study their properties 43 



52(1 . In this work, we 



focus on the productions of P-wave charm states, bottom states and bottom-strange states 
in exclusive semileptonic and nonleptonic decays. 

Since the discovery of Pso(2317) j^, the heavy-hght orbitally excited states stimulated 
continued interesting attentions. There are some special characters of these excited states, 
for example, the mass of Pso(2317) is much smaller than the prediction of the relativistic 
quark model jsl which has been very successful, and it has a narrow decay width. Though it 
is believed to be the orbitally excited state of P^ by most of the physicists now, there have 
been some arguments about its nature. It can be a conventional cs state {7,8], four-quark 
state [9|, or molecular state since it is just above the threshold of P^tt and DK 10[, etc.. 

In the family of excited heavy-light states in the conventional quark model, P^o, P^i, 
Psi and Ps2 are the orbitally excited states of P^ and they are the sb (P*j) system. We 
know little about them since only two candidates of them are observed in experiments. The 
CDF collaboration reported their observations, P^i with mass M{Bsi) = 5829.4 ± 0.7 MeV 
and Bs2 with mass M(P,2) = 5839.6 ± 0.7 MeV in 2008 [U.]- Later the DO Collaboration 
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confirmed the existence of 5^2(5840) witli mass M{Bs2) = 5839.6 ± 1.1 ± 0.7 MeV and 
indicated tliat 5^1 (5830) was not observed witli available data [l^ . 

Different from the lack of data in experiments, there are a lot of theoretical efforts to 
investigate the properties of the B*j system. For example, the mass spectroscopy had been 
estimated by the model of HQET [13|, relativistic constituent quark models [Tl, Il4j-|l7l| and 



lattice QCD [l8|. The strong decays of Bgo, B'^^, Bgi and Bs2 are also studied by many 
authors, these studies helped us not only to find another two states Sso and S^^, but also 
to estimate the full decay widths of these B*j states. 

There are large discrepancies between the existing results of the different models, which 
are shown in the section of numerical results. More careful study is needed, especially in the 
relativistic models, because the relativistic corrections are large for excited states. In this 
Letter, we will study the strong decays of BgQ, B'^^, Bgi and Bs2 by the improved Bethe- 
Salpeter(B-S) approach which is a relativistic method based on a relativistic four- dimensional 



wave equation 



20 



2l|. In this model, the B*j are bound states composed of quark s and 



anti-quark b, with an angular momentum L = 1, so they are orbitally excited states and also 
called P wave states. The quantum numbers of these P wave states are 0"^ [Bgo), 1"*" 
1+ (Bsi) and 2+ (-852), the allowed strong decay modes are 0"'' — > 0^0^, 1+ — )• 1^0~, 
2+ — > 0^0^ and 2+ — > 1^0~, while other strong decays of P wave in the final state are 
ruled out by the kinematic possible mass region. For the same reason, we have checked 
that in the final states of the allowed strong decays the pseudoscalar 0~ state must be the 
light meson {K, vr), and the other one is a heavy meson {B, B*, Bg, B*). Using the reduction 
formula, PCAC relation and low energy theorem, we got the strong decay amplitude 22| . 
which is a function of the transition matrix element between two heavy mesons. We will 
adopt this method to calculate the transition matrix element by the improved B-S method 
in this Letter. 

Similar to the cs system, the sb (or bs) system is the bound state composed of a heavy 
quark and a light quark. Since the heavy quark b is much heavier than the light quark s, the 
heavy-light mesons can be characterized by the spin of heavy quark Sg, the total angular 
momentum of light quark jg = Sg + L, and the total angular momentum J = + jg. For 
L = 0, the jg = \ H doublet, there are two states with = 0^, 1~; for L = 1, there are 
two degenerate doublets: jg = |^ 5 doublet and = ^'^ T doublet, with the corresponding 
quantum numbers = 0+, 1+ and = 1^,2+, respectively. Bgo and B'^^ are S doublet 
which are still missing in experiments, -Bsi(5830) and -Bs2(5840) are T doublet and have 
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been observed in experiments. Obviously, there are two 1"*" states: B'^^ and 5^1 (5830), we 
use the notations ^'^ and |^ to describe them respectively. 

Recently, we have resolved the instantaneous Bethe-Salpeter equation, which is also called 
Salpeter equation, and obtained numerical relativistic wave functions for different J^^'^^ 
states 13]. We also give an improved formula of the transition matrix element which 
is based on the Mandelstam formulism and the relativistic Salpeter wave functions. The 
corresponding transition matrix element is valid for any recoil momentum whenever it is 
large or small, and we have proven that this transition matrix element is gauge invariant 
when it is necessary [24^]. So in this Letter, we will use the improved B-S method to calculate 
the strong decays of the orbitally excited heavy-light states Bso, -B^i, Bgi and Bs2- According 
to the estimated masses theoretically, Bso and B'^^ have small masses, we calculate the isospin 
symmetry violating decay processes B^q — t- BsTi and 5^^^ — )■ i?*7r through 7r° — r/ mixing and 
get small widths. Considering the uncertainties of the masses, for BgQ and B'^^, we also 
calculate the strong decay channels: BgQ — t- BK and B'^^ — )■ B*K. For B^i and Bs2i as 
they have higher masses, the Okubo-Zweig-Iizuka (OZI) rule allowed decays Bsi — > B*K, 
Bs2 BK and ^ B*K are permitted, in fact B^i and 5^2 are observed through these 
decay channels. 

The Letter is organized as follows. In Sec. II, we introduce the Bethe-Salpeter equation 
and the Salpeter equation. We show the corresponding wave functions which can be obtained 
by solving the Salpeter equation in Sec. III. The method for calculating the transition matrix 
elements of corresponding decays is shown in Sec. IV; Sec. V show the formulations of the 
decay widths. Then we show our numerical results and discussions in Sec. VI. 




6 b 
FIG. 1: Strong decay of Bs2 ^ B^K^ 



4 



II. INSTANTANEOUS BETHE-SALPETER EQUATION 



In this section, we briefly review the Bethe-Salpeter equation and its instantaneous one, 
the Salpeter equation, and we introduce our notations. 



The Bethe-Salpeter (BS) equation is read as [20|: 



(^1 - mi)x{q){^2 + m2)=tj (^^(^> q)x{k) , (1) 

where x(q') is the BS wave function, V{P, k, q) is the interaction kernel between the quark 
and anti-quark, and pi,p2 are the momentum of the quark 1 and anti-quark 2. The total 
momentum P and the relative momentum q are defined as: 

mi 



Pi = aiP + q, «! 
P2 = a2P - q, 02 



nil + ^2 
m2 



nil + "^2 

We divide the relative momentum q into two parts, q\\ and q±, 

q' = qi; + ql, 

gf ^ (P ■ q/M')P^ , ^ - gf . 

Correspondingly, we have two Lorentz invariant variables: 



(ip = ^ ^ ^T = \Jql-q^ = \f-qi ■ 

When P= 0, they turn to the usual component go and |g|, respectively. 

In instantaneous approach, the kernel V{P, k, q) takes the simple form 2l| : 

V{PXq)^V{k^,q^). 

Let us introduce the notations fp{qj_) and ri{qj_) for three dimensional wave function as 
follows: 

dqp^ 



VpWl) = ^ J -^x{q^\,q±] 



(2vr) 

Then the BS equation can be rewritten as 



xiq\\,q±) = Siipi)7]iq±)S2iP2) ■ (3) 
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The propagators of the two constituents can be decomposed as: 



Si{Pi) 



+ 



J{i)(lp + CiiM — Ui + ie J{i)qp + cciM + cui — ie 



(4) 



with 



1 



2uj. 



P 



(5) 



where i = 1,2 for quark and anti-quark, respectively, and J{i) = (— 1)*+^. 
satisfy the relations: 



Here Afp{qj 



KM + KM 



r 



Introducing the notations '^p'^{q±) as: 



A 



tp\ 



KM-^Kpi 



(6) 



^tHqi.) ^ KpMj^Vp{qA.)j^Kp{ 



(7) 



and we have 



¥'p(g±) = ^p+(g±) + </'p {qi.) + Vp^{qi-) + Vp iq±)- 

Using contour integration over qp on both sides of Eq. (I3l), we obtain: 

Kp{q±)vp{qi.)Kpiq±) KMvpMKM 



^p{qi-) 



(M — ui — U2) 



{M + UI + U2) 



and the full Salpeter equation: 



{M-ui- a;2X+(g±) = A+ (gx)r/p(g±)A+ (gx) , 
{M + ui + U2)Vp' {qi_) = -Aipiq±)Vp{q±)KM ' 

= = . (8) 

For the different J^^ states, we give the general form of the wave functions (we will talk 
about them in Sec. III). Using the last two equations in Eq. (j8]), we can reduce the wave 
functions, then solve the wave functions by the first and second equations in Eq. ([8]) to 
get the wave functions and mass spectrum. We have discussed the solution of the Salpeter 
equation in detail in Refs. S] • 

The normalization condition for BS wave function is: 



27r2 



2Pn . 



(9) 
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In our model, Cornell potential, a linear scalar interaction plus a vector interaction is 
chosen as the instantaneous interaction kernel V: 

4 ct 

V{r) = Vs{r) + Vo + %® 7°K(r) = Ar + H - % ® 7%- , (10) 

6 r 

where A is the string constant and as{q) is the running coupling constant. In order to fit 
the data of heavy quarkonia, a constant Vq is often added to the scalar confining potential. 
We see that Vy{r) diverges at r = 0, in order to avoid the divergence, a factor e~"^ is added: 

K(r) = -(l-e— ) , K(r) = -^%-"^. (11) 
a 6 r 

It is easy to show that when ar <C 1, the potential becomes the original one. In the 
momentum space and the rest frame of the bound state, the potential reads: 

where the running coupling constant as{q) is chosen as: 

127r 1 



33-2Nflog{a + J^) 



With this equation and parameters shown in Sec. VI, one can find that as{mc) — 0.38 
[Nf — 3), as{mb) — 0.26 {Nf — 4), and Nf — is chosen for hq system in this Letter. The 
constants A, a, Vq and Aqcd are the parameters that characterize the potential. 



III. RELATIVISTIC WAVE FUNCTIONS 



In this section, by analyzing the parity and possible charge conjugation parity of corre- 
sponding bound states, we give the formulas of the wave functions that are in relativistic 
forms with definite parity and possible charge conjugation parity symmetry. 



A. Wave Function for ^Sq {0 ) state 

The general form for the relativistic wave function of a pseudoscalar meson with the 
quantum number — 0~ (or J^^ — 0""+ for an equal-mass system, a qq quarkonium) can 
be generally written as eight terms, which are constructed by P, qp and gamma matrices. 
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because of the instantaneous approximation, four terms with P-Qp^ become zero, the general 
form for the relativistic Salpeter wave function of a pseudoscalar state = 0^ (or J^^ 
0^^) can be written as 



V0-{q.J = [/l(g,J/' + /2(gpjM + /3(g,J +/4(g^J^^]75, (13) 

where M is the mass of the pseudoscalar meson, and fi{qp^) are functions of — . Due to 
the last two equations of Eq. (|8]): yj^r = ip^^ = 0, we have: 

MQpJ = , h{qp) = . 14 

^ 7712^1 + miUJ2 ^ m2UJi + 'miU2 

Then there are only two independent unknown wave functions /i(q'p^) and f2{qp^) in 
Eq. ([H: 

^oAQpJ = [h{qpjr + f2{qpjM-f2{qpJ4p^^^^^^^^p^ 

+h(<i.j4., p (15) 

The numerical values of radial wave functions /i, /2 and eigenvalue M can be obtained by 
solving the first two equations of Salpeter Eq. ([8]). 

One can check that in Eq. (fT3|) . which we wrote as the wave function for = 0^ (or 
J^*^ = ^) state, all the terms except the one with have positive charge conjugate 
parity, while /a term has negative charge conjugate parity. When we consider the constraint 
relations, for equal mass system, ui = U2, so = (Eq. f|T4|) ). then the whole wave function 
has positive charge conjugate parity, that is ^ state. 

In our calculation, we obtain the numerical values of wave functions in the center-of-mass 
system of the bound state, so q\\ and q± turn into the usual components (go, 0) and (0,g), 
coi = {mf + (fY^"^ and UJ2 = {m\ + cfY^"^. Then the normalization condition reads: 

— --^4/i/2M I ■ + ■ + — ■ —\ = 2M. (16) 

The numerical values of the right sides of the first two equations in Eq. ([8]) are compa- 
rable, but since M — uji — UJ2 <^ M + uji + UJ2 ioi bound state, we know that the numerical 
value of ^^^{q) is much larger than that of {q). So in the past, authors made a fur- 
ther approximation of the Salpeter equation, deleting the others in Eq. ([8]) except the first 
equation which is about the positive wave function ip'^^{q). This seems a reasonable approx- 
imation since ^^^{q) is dominant, but we point out that, we can delete the term of ip (g). 
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but that should be done after we solve the full Salpeter equation, otherwise we obtain a 
non-relativistic wave function, not a relativistic one. Since with the further approximation, 
only one equation is left, then only one unknown wave function can be solved, we have 
to choose /s = /4 = 0, and /i = /2 in Eq. f|T3|) . then the wave function Eq. f fT3|) becomes 
ifQ- = + M)'~f5, this is well known Schrodinger wave function for a pseudoscalar. So in 
our calculation, we solve the full Salpeter equations Eq. (|8]), not only the first one in Eq. ([8]). 

According to the Eq. ([7]) the relativistic positive wave function of pseudoscalar ^5*0 state 
{B or Bs in this Letter) in the center of mass system can be written as Q]: 



75, (17) 



where the bi {i = 1, 2, 3, 4) are related to the original radial wave function /j, quark mass 
rrii, quark energy Wi {i = 1, 2) and meson mass M: 

&i = TT fiio) + f2{q) : ] ,b2 = ■ ,&3 = ■ ,b4 - 



2 \ Wi + W2 J ' nil + 1712' miW2 + m2Wi' {miW2 + m2Wi) 

Inserting the expressions of (pQ^{q) in Eq. f lT7|) and corresponding ip^rlq) (which can be 
easily obtained by ip^- = ip^- — ^q-) into the first two equations of Eq. ([8]), we get two 
coupled integral equations (the explicit expressions can be found in Eqs. (24-25) in Ref. 23|). 
By solving them, we obtained the numerical values of wave functions /i, /2 and eigenvalue 
M. 



B. Wave Function for (1 ) state 



Because of the instantaneous approximation, instead of 16 terms, the general form for 
the relativistic wave function of vector state = 1^ (or J^*-^ = 1 for quarkonium) can 
be written as eight terms, which are constructed by P, q, e and gamma matrices 25l ]: 



(A) 



fi{q±) + j^f2{Q±) + j^Mi±) + ^]^f^y 



+ M/^)/5(gx) 



+ /'^ff6{q±) + {4J'^ - q± ■ e^'^)Mq^) + ^(f/'Vx " A± " s^'^)fs{q±), (18) 

where is the polarization vector of the vector meson. One should note that we use the 
same notations of the radial wave functions fi for pseudoscalar and vector mesons, but they 
are different. It should be indicated that we will use them for other states (see below), but 
we remind the readers that their numerical values are different for the different states. 
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The equations ip^- {q±) = (p^-{q_L) = give the constraints on the components of the 
wave function ipi~{q±), so we have the relations 

. . N _ + M'^f5{q±)] (mima - ujiuJ2 + qI) . , ^ _ Hq±)M{-uji + ^2) 

Ji[(l±) — 777 : r-2 ) h[(l±) — — 7 : ^ — , 

M[mi+m2)q± [miU2 + m2UJij 

ft ^ _ [■~?i/4(^^) + M'^hi.^lA-)] ijniUJ2 - m2UJi) f&{,q±)M{uJiUJ2 - mims - q\) 

Then there are only four independent wave functions f3{q±), fi{q±), f5{q±) and fe{q_L) in 
Eq. 

One can check that in Eq. f[T5|) . all the terms except those with /2 and /y are negative 
under charge conjugation operation, while the terms with /2 and fj are positive. Applying 
the constraint relations, for equal mass system, we found the terms with /2 and fi disappear, 
then the whole wave function has negative charge conjugate parity, that is 1 state. The 
similar relations hold for the following P wave states, so we will not show the details again. 

Wave functions /3(g), f^iq), f^i^q) and f&{q) will fulfill the normalization condition: 



dq IQUJ1UJ2 r M'^ U1U2 - mim2 + (f 

•j/5/6 ; r 



(27r)^ 3 [ miUJ2 + m2iOi {mi + m2){uJi + UI2) 



j. = 2M. 
(19) 



The relativistic positive wave function of '^Si state {B* or 5* in this Letter) can be written 

as: 

^+^{q) = h ^('^ + b2 f + hU± - qi. ■ £(")) + hi^ A± ■ ^^'^) 

+q± ■ e("n&5 + k^ + br^^ + bs ii. f), (20) 

where we first defined the parameters which are functions of i^Si wave functions): 

r r i^iwi + W2) ,^ ..^(mi+ms) ..^(mi + m2) 

ni = h[q) - feiq}-, ■ r,n2 = f5{q) - f&[q)-, ■ r,ri3 = h[q) + fi{q)- 



(mi + m2) ' {wi + W2) ' {wi + W2 

then we defined the parameters hi which are functions of fi and n,: 



M , M (mi + m2) , M {w2 - wi) , 1 {wi + W2) 

^1,02 = 7 -^1,03 = 7^1,04 = --/ 9T' 

2 2 [wi + W2) 2 [miW2 + m2Wi) 2 (wiW2 + mim2 - qi) 



h = —7 , rMWgln3),6e = ~ ^(MV+gln3), 

2M (W1W2 + mim2 + g^) 2M^ (wiW2 + mim2 + g^) 

, '^s Mq)M 1 W1+W2 , ,^ W1 + W2 



2M (mit(;2 + ?7i2Wi) ' 2M2 mi + m2 (mi + m2)(wiW2 + mim2 — g^^) 

Similar to the method in last subsection, where we obtained the wave functions and 
eigenvalues for pseudoscalar states, inserting ip^^{q) and corresponding ip^J{q) into the first 
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two equations of Eq. (|8]), we obtained four independent coupled integral equations (Eqs. (37- 



40) in Ref. 17|), by solving them, we obtained the numerical results of mass spectra and 



wave functions. For other states, see below, we will not show the details of how to solve 



them, interested reader can find the details elsewhere, for example, in Ref. 17| 



One should also note that we cite same notations of m^, Wi, M, P, and bi as used in last 
subsection for pseudoscalar meson, but they are different for different states. And we also 
use the same notations for other mesons, like the following P wave states. 



C. Wave function for ^Pq state 

The general form for the relativistic Salpeter wave function of ^Pq state, which = O"*" 
(or J^'-' = O^"*" for equal mass system), can be written as [26] : 

^o+(g±) = fi{q±)^+ + hiq±)M + hiq^)^- (21) 

The equations ipQ+{q±) = ipQ+{q_L) = give the constraints on the components of the wave 
functions, so we have the relations 

nf X fi{q±)ql{mi + 7712) ^ /2(g±)gi(wi - ^2) 

JsWxj = — , — \ — 2T' •'^Wj 



M{ujiuj2 + mim2 + qi) M{7niuj2 + 7n2U}i) 

Then there are only two independent wave functions /i(g±) and /2(i?±). From Eq. ([H]), we 
obtain two coupled integral equations, by solving them, we obtain the numerical results of 
mass spectra and wave functions. 

The normalization condition for the ^Pq wave function is: 

dq 16fj2U^,U2q^ ^ 2M. (22) 
The relativistic positive energy wave function of Bgo (^Pq) can be written as: 

¥'o/(^ = «i(#± + 02^^^ + % + 04^), (23) 

where the parameters are functions of /i and /2 (O"*" wave function) and are defined as: 

1/ mi + m2\ ^1 + 1^2 2 {W1 + W2) {m2Wi-7niW2) 
«i = o tiKO) + J2[q) \ ,02 = \ , % = g± \ , «4 - 



2 \ 'Wi + W2 J mi + m2 miW2 + 7n2Wi (mi + 77^2) 
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D. Wave function for ^Pi state 



The general form for the Salpeter wave function of ^Pi state, which = l"*" (or J^*^ 



1++ for equal mass system), can be written as 



26|: 



where e^'^-' is the polarization vector of the ^Pi state. 
The constraint equations give us the relations: 



fi{q±)M{miU2 - m20Ji) f2{q±)M{-uiu:2 + mim2 + q\) 
— ; U\q±) — 



q]_{uji + UJ2) ' gi(mi + m2) 

The normalization condition for the ^Pi wave function is: 

r dq 32fif2UiU2{^^i^^2-mim2 + ^) _ 
J (27r)3 3(mi + m2)(u;i + ^2) 

The relativistic positive energy wave function of ^Pi state can be written as: 



(25) 



^t^^iq) = te^^^^P'qle^^^^ailMY + 027^ f + 4i. + a^l^ F (26) 
where the parameters aj are functions of /i and /2 {^Pi wave function) and are defined as: 



2 V 7711 + 1712, 



mi + 7n2 M{wi — W2) 

• ^3 = : ) ct4 



Wi + W2 



77I1W2 + 77I2W1 ' 



(mi + m2) 
miW2 + m2Wi 



E. Wave function for ^Pi state 

The general form for the Salpeter wave function of ^Pi state, which = 1^ (or J^*" 
l"*"" for equal mass system), can be written as 26| : 



V5iPi(g±) =q±-e 



(A) 



+ f2iq±)^ + h{qi)^ + /4(9±)0 



75, 



where is the polarization vector of the ^Pi state. 
The constraint equations provide us the relations: 

/i(g±)M(mi - m2) 



-, f4{q±) 



f2{q±)M{ui+L02] 



(^1^2 + mim2 - gi) ' •'"^'^-^^ (mia;2 + "^2^1) 

The normalization condition for the ^Pi wave function is: 

dq 16/1/2^1^2?^ 



(27r)3 3(miCJ2 + m2a;i) 



2M. 



(27) 



(28) 
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The relativistic positive energy wave function of ^Pi can be written as: 



75, 



(29) 



where the parameters are functions of /i and /2 (^Pi wave function) and are defined as: 

1^1 — W2 mi + m2 



ai = o /i^ + /2 , , 

/V mi+ 1712 



mi + 777.2 

. a2 = ; , as 

Wi + 7772 



m27yi + mi7i'2 

1+ „v.^ 3 



, a4 



7r727i'l + mi7y2 



The wave functions of two physical 1^ states (or and 1"^) are the mixing of ^tp^{q) 
and iptp^{q), see Eq. (15^ below. 



F. Wave function for P2 state 



The general form for the relativistic wave function of tensor = 2^ state (or J^'" = 2'^'^ 
for equal mass system) can be written as 



(30) 



where e''^} is the polarization tensor of the 2+ state. The constraint equations give further 

[qifM + M^um (^1 + ^2) - M^um^i - ^2) 



relations: 



f2{q) 



M{miU2 + m2U)i) 
[qlfM-M'Mq)] {001-002) 



M{miUJ2 + 7772^1) 



friq) 



f,{q)M{coi-0O2] 



fs{q) 



h{q)M{ooi + 002) 



(31) 



7r7ia;2 + m200i miU2 + 7772^1 

Only four independent wave functions /3(g), f-i^q), h{.q) and /6(i?), the numerical values and 
the bound state mass M can be obtained by solving the full Salpeter equation. 
These four independent wave functions fulfil the normalization condition: 



16 001IJO2 q 



7fi 



+/4 /5 q 



{2tc)^ 15{miL02 + m200i 

{rrii + 7/72) (7/72^1 — 7r7iu;2 



5 + 



(7771 + 7772)(7772Wl - 7771^2) 

u;iu;2(wi + U2) 



2 + 



2g73 /4^ + /6 



2M. (32) 



iOiC02{cOl + OO2) 

The relativistic positive energy wave function of Bs2 (^-^2) can be written as: 
y^sV^ = + a2— + 03— + «4^^^] + 7'^[a5 + + 07— + ag^^jl^S) 



13 



similar to 1 state, we first defined nt as: 

ni = /3(g) + U{q) : , n2 = h{q) - feiq) : , 

Wi +W2 mi + 1712 

then we defined the parameters af. 

^ ^ (wi + W2)ql ^ ^ if5iq)w2 - fe{q)m2)M 
^ 2M{miW2 + m2Wi) ^ 7711^2 + ^2^1 ' 

(7771 - m2)ql {f6iq)W2 - h{q)m2)M 

0-2 — — TTT T*^! ~r , 

2M (77717^2 + 'T727i'l) 77717772 + 

1 , /6(glM^ 1 77^1+7772. ^ fMM^ 



2 77717772 + 7772 7771 2 TTii + 7ri2 7rii77'2 + m2Wi 

M M(mi+m2) M^(wi-W2) M^(mi+m2) 

as = -7rn2, ag = — ■ -r72, aj = — ■ rr72, ag = — ■ rr72. 

2 2[wi + W2) 2{m2Wi + -1711102) 2{m2Wi + 771111)2) 

IV. TRANSITION MATRIX ELEMENT 

In this section, we show the method to formulate the transition matrix element, which is 
general for all the decay channels in this Letter. 

Considering the limitations of phase spaces, there are seven dominant strong decay chan- 
nels for Blj states: B^q 5°7r°, B^o ^ BK, B'^-^ Bfn°, B'^^ B*K, B^i B*K, 
Bs2 BK and 5,2 ^ B*K (where BK = B+K' + and B*K = B*+K- + B*^K^). 

Since all the light mesons in final states are pseudoscalars, we can give a unique formulation 
of the transition matrix element for these seven decay channels. By using the reduction 
formula, PCAC relation and low energy theorem, taking the channel Bs2 BK as an 



example, see Figure 1, the corresponding transition matrix element can be written as 22l |: 



T{Bs2 ^ = — < B\Pf)\d^,^,s\B,2{P) >, (34) 

/ K 

where fj^ is the decay constant of pseudoscalar K meson, P' is the momentum of K. The 
contribution of the light pseudoscalar is reduced to a factor -7-^, then the main part of the 
calculation in Eq. is to calculate the transition element < i?°(P/)|(i7^75s|i?s2(-P) >■ 
If we further choose the instantaneous approach, according to the Mandelstam formalism 



281], at the leading order, the transition matrix element can 



)e written as an integral equation 



of the corresponding initial and final state wave functions 



{B%Pf)\dj,j,s\B,2iP)) = I (||^2^r{^2o-^(g07.75^^:(g1^} , (35) 
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where P and M are the momentum and mass of initial state Bs2', Q and q' = q + ^^^^ x Pf 
are the relative momenta of quark and anti-quark in the initial state 3^2 and the final state 
respectively, which are defined as q = ps = —pb and q' = ^^^^ ^ Pf — Pb in the center 
of mass system of initial state Bs2] V'^^(Q'-l) ^^id <^'^^(5'j_) are the positive energy wave 
functions of Bs2 and 5°, which are given in last section. 

In our model, improved B-S method, which is based on the constituent quark model, 
we give the forms of wave functions by considering the quantum number or J^^ for 
different states, and these states in our model are labelled as ^Sq (0~ state), ^Si (1~ state), 
^Pj (J = 1, 2, 3) (0+, 1+, 2+) and ^Pi (1+). For the unequal mass system, the ^Pi and ^Pi 
states are not physical states, the two physical states |^ and I"*", which are the mixtures of 



them, can be expressed as [l3|, |29|-|31|: 







1 + 






>= 


'2 


> 


\Bsi 




3 + 




>= 


'2 


> 









:ose\^Pi> +smefPi>, (36) 

where 6 is the mixing angle and 6 ~ 35.3° in the heavy quark limit. 

The strong decay amplitudes can be described by the strong coupling constants, they are 
defined as: 

T{BsO — > BsTl) = GbsoBsTt, 

T(5,o -> BK) = Gs^^BK, 

TiB',, ^ B» = Gs'^.Ba^?^ ■ ^)(^^'^ • + FBi^BtM''^ ■ 
T{B'^, ^ B*K) = Gb'^^b^M''^ ■ v)ie^'^ ■ v') + " e^^^). 

T{B,, ^ B*K) = GB^^B^Rie?^ ■ v)is^'^ ■ v') + " ^^'^), 

T{B,2 ^ BK) = GB^,BK^^i}v'^v'\ 

T{B,2 ^ B*K) = zGB^.B^Rsi'Jv'^e^'^"'''^'^'^, (37) 

where Gb^qBstt^ ■ ■ ■ Gb^2B*k are the strong coupling constants; v = ji and v' = jfj are four- 
velocities of initial state and final state; e is the polarization vector of initial state B'^^ or 
Bgi, 61 is the polarization vector of final state B* or B*, e^^, is the polarization vector of 
initial state Bs2', P, Pf are the four- momenta of initial and final heavy states, respectively. 
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V. DECAY WIDTHS 



The decay channels of S,o ^ BK, B'^^ B*K, B,i B*K, B,2 ^ BK and ^ B*K 
are OZI rule allowed, through the Eq. fl34|) and Eq. f p5|) the calculations of corresponding 
decay widths are straightforward: 



71 



BsoBK 



Bs2B*K 



IP 



/I 



247rM2 



^ |t(5,i ^ B*K) 



\Pf 



A, A' 



Bs2BK 



A,A' 



71 



407rM2 



IP 



/I 



407rM2 



^ |t(P,2 ^ p*ir) 



A, A' 



(38) 



The decays of BgQ — )• PjTt and — )■ P*7r violate the isospin symmetry which only occur 
through TT^ — 7] mixing p^]. According to Dashen' theorem 32|, the decay widths can be 
written as: 



IP, 



BsoBsTT 



f\ 



87rM2 



ml — ml 



IP 



/I 



247rM2 



E 

A,A' 



T(P;i ^ P:7r)t., 



mi — 



(39) 



where t^r.^ =< 7r''|'H|r7 > is the tt^ — rj transition matrix, m,^ and are the masses of vr and 
77. The chosen value of t-j^ri = —0.003 GeV^ 32(| is very small, which result in narrow decay 
widths of these two channels. 



VI. NUMERICAL RESULTS AND DISCUSSIONS 

In order to obtain the masses and wave functions of initial and final states, which are used 
to calculate the transition matrix elements and decay widths, we solve the instantaneous 



20I, 



B-S equation 
Aqcd = 0.27 GeV, a 



2l| . and the parameters are chosen as: a = e = 2.7183, A = 0.21 GeV^, 
0.06 GeV, mb = 4.96 GeV, =_0.5_GeV, jn„ = 0.305 GeV, 



HQ 
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-|27l|. The numerical 



ma = 0.311 GeV and Vq for the B-S kernel as same as Refs. 
values of these parameters are obtained by fitting data, for examples: Mb^^ = 5.830 GeV, 
Mb^2 = 5.840 GeV, Mb = 5.279 GeV, Mb- = 5.325 GeV, Mb, = 5.366 GeV, M^* = 5.415 
GeV [331. 



With this set of parameters, and varying all the input parameters simultaneously within 
±5% of the central values, we obtain the light doublet masses: 



Mb,^ = 5.723 ± 0.280 GeV, M^/^ = 5.774 ± 0.330 GeV. 



(40) 
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TABLE I: The strong decay widths of BgQ and B'^i in units of MeV. Here we only calculated the 
cases that the upper limit values of masses of Bso and B'^^ are chosen, and we take 6 = 35.3°. 



Mode 


Ours 


B,o ^ B+K- 
Bso ^ 


67.1 
67.0 


B^i ^ B*+K- 


37.7 
37.6 



The other input parameters in this Letter are the masses of hght mesons: M^r = 0.139 GeV, 
Mj = 0.548 GeV, Mk = 0.494 GeV, and the decay constants = 130 MeV, fx = 156 MeV 



33|. 



With these parameters, we solve the instantaneous Salpeter equations for different states, 
and obtain the masses and relativistic wave functions. We calculate the transition matrix 
elements by the wave functions numerically, and get the strong coupling constants: 

GbsoBs. = 17.5 ± 2.5 GeV; = -16.3 ± 2.2 GeV, F^.^bj. = 17.2 ± 3.3 GeV; 

GbsiB*+k~ = -2.39 ± 0.36 TeV, Fb,,b'+k~ = 0.320 ± 0.065 GeV; 
Gb^,b*oko = -2.48 ± 0.37 TeV, Fb^^b*oko = 0.227 ± 0.050 GeV; 
Gb,2B+k- = 1-96 ± 0.26 TeV, Gb,,bo j?o = 1.97 ± 0.26 TeV; 

Gb,2B*+k- = 2.33 ± 0.65 TeV, Gb^2b*oko = 2.35 ± 0.65 TeV. (41) 

Considering the uncertainties of the masses, for mesons BgQ and B'^^, the upper limits of 
the masses are Mn ~ 6.00 GeV and Mn' ~ 6.10 GeV, which are above the threshold of 
BK, B*K, so with these upper limit values of masses, we also calculated the OZI allowed 
channels: 

Gb^oBORo = 28.6 GeV, Gb^^b+k- = 28.7 GeV; 
Gb' b'ORo = 2.34 TeV, Fb' b-oro = 6.34 GeV; 

si si 

Gb'^,B'+k- = 2.32 TeV, Fb'^^b*+k- = 6.34 GeV. (42) 

The corresponding decay widths are shown in Table L 

In Table. I, II, III and IV, we show numerical results of strong decay widths, which are 
predicted by us and other authors. For Bgo and B'^^, our results are consistent with the ones 
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TABLE II: The strong decay widths of Bso and B'^i in units of keV, where the decays of Bgo — )■ BK 
and B'g^ — t- B*K are calculated with the upper limit masses of BgQ and B'^^. And the mixing angle 
6 = 35.3° is used. 



Mode 


Ours 


[43j 


[44,45] 


[46] 


m 


[50] 


Bso BgiT 


13.6 ±5.6 


6.8~30.7 


1.54 [44] 


55.2~89.9 


21.5 




Bso ^ BK 


134 MeV 










227 MeV 


B'si ^ i?> 


13.8 ±3.6 


5.7~20.7 


10.36 [45] 


57.0~94.0 


21.5 




^ B*K 


75.3 MeV 










149 MeV 



TABLE III: The strong decay widths of Bgi and Bs2 in units of MeV. Here we take 9 = 35.3°. 



Mode 


Ours 


Bsi B*+K~ 
Bsi B*°K° 


0.028 ± 0.0075 
0.013 ±0.0036 


B,2 ^ B+K- 
Bs2 ^ B^^ 


0.83 ±0.23 
0.72 ±0.20 


Bs2 ^ B*+K- 
Bs2 B*^K^ 


0.091 ±0.052 
0.057 ±0.032 



presented by references |43j and 45] , close to the results of reference 47|], but much larger 



than the ones of 



44[, and much smaller than the predictions of 46|. Though the predicted 



masses of the P-wave -B*j states are similar for different models, the predicted decay widths 
are much different. The situation is similar in other channels. For example, our prediction 
of Bsi —7- B*K is close to the result of reference 48|], and we get the narrowest decay width. 

In our calculation, we choose a small value of mixing parameter, = —0.003 GeV^, 
which suppresses the corresponding decay widths much heavily, the decay widths are very 
small. Considering the uncertainties of the masses, the decay widths of the OZI allowed 
decay channels: Bsq — ?■ BK , B'^^ — t- B*K are larger than the results of Bsi and Bs2- This is 
because that the S doublet states (0^, 1^) decay through S wave, they usually have broader 
decay widths, while the T doublet states (1^, 2+) decay through D wave, they usually have 
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TABLE IV: The strong decay widths of B^i and i?s2 in units of MeV. Here we take 9 = 35.3°. 



Mode 


Ours 


[15] 


[48J 


[49] 


[50] 


[51J 


[52] 


B,i B*K 


0.041 ±0.011 




0.098 


3.5 


0.4 ~ 1 


0.28 


< 1 


Bs2 BK 


1.55 ±0.43 


2.6(1.9) 


4.6 




2 


7±3 


1 


Bs2 ^ B*K 


0.148 ±0.084 


0.07(0.05) 


0.4 


3.2 


0.12 




< 1 



narrower widths, though the decay channels are OZI rule allowed ones. 

The large discrepancies between the results of different models may be caused by the 
small phase space of transition channels, and the results are very sensitive to the masses of 
corresponding mesons. We take into account the errors by varying all the input parameters 
simultaneously within ±5% of the central values. One can see that in Table. II, we get 
relatively large errors even we change the parameters in such small regions. In Table. Ill 
and IV, the phase spaces of the decays of -B^i and are too narrow, the decay widths 
depend heavily upon the masses of initial mesons. If we change the masses of initial mesons, 
there will be large errors, even one order larger than the original value. So we change all the 
input parameters except the masses of initial mesons which are fixed to the experimental 
data to calculate the errors for B^i and Bs2- 

In conclusion, through the improved B-S method, we predict the masses of the orbitally 
excited states BgQ and B'^^, calculate the strong coupling constants Gb^oBs-k, Gb^qBKj Gb'^^b*-k-i 
Fe'^^bi-k, Gb'^^b*k, Fb'^^b*k, iB*K^ FbsiB*k^ Gbs2BK^ Gbs2B*k^ aiid obtain the strong 
decay widths of B,o ^ B.ir, Bso BK, B'^^ ^ B*7r, B'^^ B*K, B,i B*K, Bs2 ^ BK 
and Bs2 — ?• B*K, which are useful to find the unobserved states and to estimate the full 
decay widths of these orbitally excited states. 
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